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We calculate the third coefficient of the lattice beta function in QCD with Wilson fermions, extending the pure 
gauge results of Liischer and Weisz; we show how this coefficient modifies the scaling function on the lattice. We 
also calculate the three-loop average plaquette in the presence of Wilson fermions. This allows us to compute the 
lattice scaling function both in the standard and energy schemes. 



1. INTRODUCTION 

In lattice formulations of field theories the rel- 
evant region for continuum physics is the one 
where scaling is verified. Scaling in the continuum 
limit requires that all RG-invariant dimension- 
less ratios of physical quantities approach their 
continuum value with non-universal corrections 
which get depressed by integer powers of the in- 
verse correlation length. Tests of asymptotic scal- 
ing involve the relationship between the bare cou- 
pling and the cutoff as an essential ingredient. 
This relationship is expressed in the bare lattice 
P function which reads, in standard notation: 



P (3o)=-a— 
da 



= -&0.90- 



■b{gl-b^gl-...{l) 



The first two coefficients in Eq. (1) are scheme- 
independent and well known. Our purpose is to 
compute the third coefficient, b\ , in SU (N) Yang- 
Mills theory with Nf species of Wilson fermions. 

M. Liischer and P. Weisz [Q have calculated 
this coefficient in the pure gluonic case, using a 
coordinate space method for evaluating the lat- 
tice integrals. We have verified their result [Q, 
using a different technique where firstly the inte- 
grands are Taylor expanded in powers of the ex- 
ternal momenta and then computed term by term 
with the introduction of an IR regulator. In the 



present work, we are extending our calculation to 
include contributions from Wilson fermions. 

We also present here the calculation of the 3- 
loop coefficient u>3 of the average plaquette, in 
QCD with Wilson fermions. In the absence of 
fermions, this quantity was calculated by some of 
the present authors H. 

Our results for b% and W3 allow us to study the 
corrections to asymptotic scaling in Monte Carlo 
simulations of dimensionful quantities, for both 
the standard and energy schemes. 

The involved algebra of lattice perturbation 
theory was carried out by making use of an ex- 
tensive computer code developed by us in recent 
years Q ; for the purposes of the present calcula- 
tion, the code was extended to include fermions 
and form factors. 

2. CALCULATION OF 

We performed this calculation using the back- 
ground field method. In order facilitate compari- 
son we have adopted the notation of Ref . jl| . 

We set out to compute d,2{^a) in: 

"MsO") = a n + di{[ia)al + d 2 {li.a)al + . . . , 



= .9o/ 47r > "msM = 9 2 <A>/ 47r 



(2) 



*Presented the paper. 



As we shall see, b\ will follow immediately from 
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this quantity. d2{pa) can be expressed as: 
d 2 {pa)={A^)\{v R X \p)-v^{p)) 2 -vf{p)+V^{p) 

- X(du^ (p)/d\) (p) - w £) (p))} A=Ao (3) 

in terms of the gauge parameter A, and of w(p) 
and i/(p) (the gauge- and background- two point 
lattice functions, respectively): 

E r(2A0) fe -p)&= - - "(pM 

E r( °' 2 ' 0) (P, -P)^= - [3(1 - u;(p)) + A ] 

oo 

Kp) = £s?" (,) (p) 

2=1 

oo 

W ( F )=E^ W (P) (4) 
i=i 

The subscript i? refers to the analogous MS renor- 
malized quantities. 

The 1-loop quantitites appearing above, 
fW.fjJ'.wfVii', are known. For the MS 2- 

(2) 

loop function v R (p) there are four new diagrams 
with respect to the pure gluonic case. We find: 

4 2) (A=1, J3)(16tt 2 ) 2 = A 2 [8p + 577/18 - 6C(3)]+ 

Nf[(-3p - 401/36) A +{p + 55/12 - 4C(3))/A] 

with: p = \n(p 2 /p 2 ). (5) 

The only remaining quantity to calculate (and 
by far the most complicated) is i/ 2 )(p). The in- 
clusion of fermions brings in 20 additional dia- 
grams (of these, two correspond to fermion mass 
rcnormalization) . 

The evaluation of these diagrams requires 
tremendous analytical effort. We have developed 
an extensive computer code for analytic compu- 
tations in lattice perturbation theory Q. In a 
nutshell, this computer code performs: Contrac- 
tion among the appropriate vertices; simplifica- 
tion of color/Dirac matrices; use of trigonome- 
try and momentum symmetries for reduction to 
a more compact, canonical form; treatment of 



(sub)divergences and extraction of logarithms by 
opportune additions/subtractions (in the case at 
hand, this step gives rise to hundreds of types of 
expressions, each containing typically hundreds of 
terms); automatic generation of highly optimized 
Fortran code for the loop integration of each type 
of expression. 

The integrals are then performed numerically 
on finite lattices. Our programs perform extrap- 
olations of each expression to a broad spectrum of 
functional forms of the type: J2i j e^lni)- 7 /D, 
analyze the quality of each extrapolation using a 
variety of criteria and assign statistical weights 
to them, and finally produce a quite reliable esti- 
mate of the systematic error. 

Several consistency checks can be performed 
on the separate contributions of each diagram: 
O(p ) parts obey several relations, and must sum 
up to zero by gauge invariance; Lorentz non- 
invariant terms ^Z^P^/p 2 must cancel; the single 
and double logarithms are related to known re- 
sults. We have performed all these checks both 
analytically and numerically; this was also a use- 
ful verification of the correctness of our error es- 
timates. 

Our final result for v^> (p) is (we take r = 1 for 
the Wilson parameter): 

v^{p) = ^ 2 \p)\ Nf=0 (6) 

+N f [(-1/N + 3N) ln(p 2 )/(167r 2 ) 2 
-0.000706(3) /A + 0.000544(5) * A] 

Also, we have: 

u^(p)-v^(N f =0,p) = J 1 \p)~J 1) (N f ^0,p) 
= N f * [ln(p 2 )/247r 2 - 0.01373219(1)] (7) 

A table of partial results will appear in a forth- 
coming publication g, where we also expect to 
provide better accuracy. 

From the relations given above, we can now 
easily construct di(pa) and e?2(/ia). Writing: 
d\{pa) = dio+dn ln(pa), d 2 (pa) — {di{pa)) 2 + 
d 2 o + $21 hi(/ia), the coefficient b\ follows imme- 
diately: 

b% = bf* + {d 20 d n - d 21 d 10 )/ (128tt 3 ) (8) 
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where 

bW = ^— - r [-57147V 3 + 7V?(-224iV + 66/iV) 

1 108(4tt) 6L > y 1 ' 

+ N f (M18N 2 - 561 - 27/N 2 )} (9) 

In particular: 

62 (iV=3, AT/=3, r=l) = -0.002284(3) (10) 

(cf. &£(jV"=3,JV/=0) = -0.00159983) (11) 

The deviation, q = (b\ - b%b )/(2bl), in the 
asymptotic scaling formula: 

aA L = cxp(-l/2fe go 2 )(&o.9o 2 )" bl/2bo [l+<ZSo + • ■ •] 
is now found to be: 

q(N=3, N f =3, r=l) = 0.3694(4) (12) 

3. THE 3-LOOP MEAN PLAQUETTE 

We have computed the free energy to 3 loops 

- In Z/V = const. + F 2 g 2 + F 3 g% + ... (13) 

in QCD with Wilson fermions. A total of 24 dia- 
grams involving fermions contribute to F3. 
The equality: 

(1 - D/N) = -{l/6V)(d\n Z/d(3) (14) 

(valid by virtue of (Sfermion) = beyond tree 
level) relates F2 and F3 directly to the quanti- 
ties W2 and W3 in the expansion of the average 
plaqucttc: 

(1 - n/N) = Wl gl + W29t + wsSo + ■ • • (15) 

We report below the final values for w 2 and W3 
at 7V=3, Nf=3, r=l. A breakdown of the results, 
along with the explicit N- and ^-dependence 
and improved error estimates will be presented in 
a forthcoming publication . 

From the plaquette coefficients, we may now 
construct the coupling constant in the "energy" 
scheme: 

9e = 9o + (w2/wi)5o + { w ^/wi)gt + ■■■ (16) 

In terms of this improved definition of the cou- 
pling constant, the j3 function reads: 

&f = b 2 + {bow 3 - biw 2 - b wl/wi)/wi (17) 



Now, the scaling violation parameter q — (b 2 — 
6f6o)/(2&o) becomes (at N = 3, N f = 3,r = 1): 

q(k = 0.1675) = 0.2103(4) 

q(k = 0.1560) = 0.1764(4) (18) 

This indeed represents an improvement over the 
standard case. We are presently working on im- 
proving the accuracy of our results, for the longer 
write-ups |5|||. 



k 


10 3 X W 2 


10 3 x W3 


0.156 


18.2522(2) 


9.15(1) 


0.1575 


17.7713(2) 


8.91(1) 


0.16 


16.9767(1) 


8.51(1) 


0.164 


15.7243(2) 


7.89(1) 


0.1675 


14.6480(6) 


7.35(1) 



Table 1 

Values of w 2 x 10 3 and w 3 x 10 3 . 7V=3, iV/=3, 
r=l, k: hopping parameter. 
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